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ABSTRACT 

We revise a magnetic buoyancy model that has recently been proposed as a 
mechanism for extra mixing in the radiative zones of low-mass red giants. The 
most important revision is our accounting of the heat exchange between rising 
magnetic flux rings and their surrounding medium. This increases the buoyant 
rising time by five orders of magnitude, therefore the number of magnetic flux 
rings participating in the mixing has to be increased correspondingly. On the 
other hand, our revised model takes advantage of the fact that the mean molecular 
weight of the rings formed in the vicinity of the hydrogen burning shell has been 
reduced by ^He burning. This increases their thermohaline buoyancy (hence, 
decreases the total ring number) considerably, making it equivalent to the pure 
magnetic buoyancy produced by a frozen-in toroidal field with B^p ^ 10 MG. We 
emphasize that some toroidal field is still needed for the rings to remain cohesive 
while rising. Besides, this field prevents the horizontal turbulent diffusion from 
eroding the contrast between the rings and their surrounding medium. We 
propose that the necessary toroidal magnetic field is generated by differential 
rotation of the radiative zone, that stretches a pre-existing poloidal field around 
the rotation axis, and that magnetic flux rings are formed as a result of its 
buoyancy-related instability. 

Subject headings: stars: abundances — stars: evolution — stars: interiors 



1. Introduction 



During their flrst ascent on the red giant branch (RGB), a majority of low-mass stars 
(those with M < 2 Mq) experience extra mixing in their radiative zones separating the H 
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burn i ng shell from the bottom of convect ive envelope (jSweigart &: Mengellll979l : ICharbonnel fc Do Nascimei 
19981 : iDenissenkov &: VandenBere]|2003l ). Despite 30 years of effort, however, the underlying 
physical mechanism is still not understood. Observationally, the RGB extra mixing mani- 
fests itself through changes of the surface abundance s of Li, C, N, and of the isotopic ratio 
i2q^ i3q correlating with an increasing luminosity (IGratton et al.l l2000l : ISmith fc Martell 
20031 ) ■ These changes are produced by the joint operation of thermonuclear reactions that 
take place in the vicinity of the H shell and a nonconvective mixing process that transports 
reaction products through the radiative zone to the convective envelope. Observations sup- 
port the idea that this mixing process starts (or, gets much more efficient) when an RGB star 
reaches a luminosity at which the differential luminosity function for a population of stars 
having the same age and chemical composition shows a prominent bump (a local pile-up of 
stars). The luminosity bump results from a temporary slowing down of the star's evolution 
caused by its structural readjustment. This happens when the H shell crosses and erases a 
discontinuity in the H-abundance profile left behind by the bottom of convective envelope at 
the end of the first dredge-up. During the first dredge-up, that occurs on the subgiant branch 
and on the lower RGB, the convective envelope grows in mass, which causes its bottom to 
penetrate the layers whose chemical composition had been altered yet on the main sequence 
(MS). This produces changes of the surface abundances of Li, C, N, and of the ^^C/^^C ratio 
similar to but by far less substantial than those incurred from the subsequent operation of 
the RGB extra mixing. 

Until recently, it has been thought that the only reason why the RGB extra mixing does 

not manifest itself below the bump luminosity is a strong gr adient of the mean molecular 

weigh t /i caused by the onset of a deep convective envelope (e.g.. lCharbonnel. Brown, fc Wallerstein 
19981 ). Any mixing mechanism has to overcome the stable thermal stratification of the radia- 
tive zone, and in the presence of a positive such mixing is correspondingly more difficult. 
In an RGB star above the bump luminosity, the H shell has already crossed the H-profile 
discontinuity, therefore the radiative zone is now chemically uniform everywhere except in a 
very close neighborhood of the H shell. This circumstance was repeatedly emphasized in the 
past. In particular, it has been used to model the RGB extra mixing with rotation-driven 
meridional circulation and turbulent diffusion. It was not until recently that it has become 
clear that rotational mixin g fails to explain the chemical elemeiit transport in the radiativ e 
zones of upper RGB stars (jChaname. Pinsonneault. fc Terndrupll2005l : iPalacios et al.ll2006l ). 
In short, this failure is due to the following main causes: firstly, rotation period measure- 
ments for young cluster stars and helioseismic data indicate that MS stars with M < 1 M© 
loose a great amount of their initial angular momentum via magnetized stellar winds and 
that they most likely become slow and nearly solid-body rotators before leaving the MS; 
secondly, the chemical element transport by meridional circulation is strongly hindered by 
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rotation-induced horizontal turbulence in stellar radiative zones (jChaboyer fc Zalinlll992l ): 
thirdly, the vertical turbulent diffusion powered by differential rotation in the radiative zones 
of RGB stars operates at a low level too because it also redistributes the angular momentum, 
thus reducing the degree of differential rotation in a self-regulating way. 



A new class of RGB extra mixing models has emerged since lEggleton. Dearborn, fc Lattanzio 



( I2OO6I ) noticed that a tiny /i-gradient inversion (V^ ~ —10 becomes visible at the outer 



tail of the H burning shell precisely at the moment when the H shell erases the H-profile 
discontinuity. This inversion is produced by the reaction ^He(^He, 2p)^He that locally re- 
duces the mean molecular weight by A/x ^ fi^AX^/G, where X3 is the ^He mass fraction. 
The mechanism can be effective in the low-mass RGB stars because their MS progenitors 
synthesize large amounts of '^He in their outer radiative cores through non-equilibrium pp 
burning. Even though this ^He-rich material gets diluted in the convective envelope during 
the first dredge-up, the radiative zone of a low-mass RGB star above the bump luminosity 
can still have X3 increased up to a value of 2 x 10~^ (the solar initial ^He abundance is 
3 X 10~^). For this mass fraction, the '^He burning leads to A/i ^ — 10~^, assuming that 
fi ^ 0.6 and AX^ f» —^3- Below the bump luminosity, the /x-gradient inversion is overridden 
by the strong positive /i-gradient built up on the MS. It shows up and may come into play 
only when the ^He burning shell, advancing in mass in front of the major H shell, finds itself 
in the chemically homogeneous part of the radi ative zone. This happens at the bump lumi- 
nosity. lEggleton. Dearborn, fc Lattanzio! (120061 ) found that a rapid mixing process occurred 
in t heir 3D simulations above this point , although the underlying cause was not identified 



see 



Denissenkov &: Pinsonneaultl l2008bl ) 



Inspired by this work, ICharbonnel fc ZahnI (j2007al ) have proposed that the /i-gradient 
inversion maintained by the ^He burning drives thermohaline convection in the radiative 
zones of low-mass RGB stars above the bump luminosity and that this is the long-sought 
physical mechanism for the RGB extra mixing. Thermohaline c onvection is a mixing process 
triggered by a double diffusive instability (e.g.. 



Vauclaidl2004[ ) 



Consider a stratified ideal 
gas with a stable temperature gradient (V = d\nT/d\nP < Vad, where "ad" stands for 
adiabatic changes) but with an unstable composition gradient (V^ < 0). If we isolate a gas 
blob and shift it up in the vertical direction then its further motion will depend on how 
fast the blob exchanges heat and composition with its surrounding medium horizontally. 
Indeed, the relative difference in density between the surrounding medium and the blob is 
Ap/p ^ Ap/p — AT/T, assuming that AP = 0. For the blob to continue rising, we need 
Ap > 0. Our assumptions about the gradients mean that Ap > and AT > in the absence 
of both heat and molecular diffusion. Because these differences grow when the blob rises, 
Ap may stay positive or it may ultimately become negative depending on the ratio = 
|V^|/(Vad — V). In our particular case, <C 1. Therefore, our idealized impermeable and 
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adiabatic blob will rise a short distance and then stop, when the accumulated difference in T 
compensates that in fi. In reality, the heat exchange, whose rate is specified by the radiative 
diffusivity K, constantly works to reduce the difference in T. On the other hand, molecular 
diffusion u^oi tries to smooth out the difference in /i. The double diffusive instability may 
therefore develop only if i^T ^ z/moi- In this case, the blob's rising speed can be estimated 
as f ~ l/rth, where / is the mean path that the blob travels before it gets dissolved, while 
''"th ~ (P/K is the characteristic thermal time scale for a spherical blob of the diameter 
d. An approximate expressi on for the thermohal i ne diffu sion coefficient can be obtained as 
Krn{l/d)'^. ICharbonnel &: ZahnI (j2007al ) and iDenissenkov &: Pinsonneault 



D 



Ivr, 



(l2008bl ) have demonstrated that the observed RGB mixing patterns can be explained by 
stellar evolutionary models with the ^He-driven thermohaline convection only ifl/d > 10-30. 
A similarly la rge par a meter ratio for thermohaline convection in stellar r adiative zones was 
postu lated by lUlrichI (Il972l ). as opposed to a ratio l/d ^ 1 advocated by iKippenhahn et al. 



(Il980l ). Besides, the double- diffusive instability has been shown to result in formation of 
elongated (large / to d ratios) structures known as "salt fingers" in labor atory experi ments 
with the saltier and warmer water overlying the fresher and colder water (jSterrull960l ). 



However, there appears to exist observational and theoretical arguments challenging this 
model. First of all, a large number of old metal-poor MS stars with M < 0.9 M©, both in 
globular clusters and in the halo field, that had accreted He- and C-rich high-/^ material from 
their evolved cluster or bi nary companions do not seena to have been thoroughly mixed by 



thermohaline co nvection (INewsham &: Terndrupl 120071 : IDenissenkov fc Pinsonneault 



2008a 



Aoki et al. I l2008h . as it would be expected even if the less efficient prescription by lKippenhahn et al.l 
Jl980h were used for -Dthc- Second, thermohaline convection is expected to be suppressed by 



the rot at ion- induced horizontal turbulence that works together with th e molecular diffusion 

to redu ce the /x contrast between the rising gas blob and its surroundings (IDenissenkov fc Pinsonneault 
2008bj). Third, strong differential rotation is predicted to hinder thermohaline convection as 
well, because the "salt fingers" may be tilted by the r otational shear so rapidly that they will 
get damped before they produce significant mixing (jCanutdll999l . and references therein). 
We anticipate that a similar effect is also produced by the Coriolis force in a uniformly rotat- 
ing radiative zone. Contrary to these expectations, a much larger fraction of Li-rich objects 
has been found among rapidly rotating (f sin i > 8kms~^ ) K giants than am ong their more 
common slowly rotating (f sin i < lkms~^) counterparts (IDrake et al.ll2002l). The Li- rich K 
giant s are low-mass stars located above the bump luminosity (ICharbonnel fc Balachandran 
2000l ) in wh ich large amounts of Li are thought to be synthesized via the ^Be-transport 
mechanism (jCameron fc Fowled 1 19711 ). To be efficient, this mechanism needs a 10 to 100 
times faster mixing t han that required to reproduc e the abundance patterns in the majority 
of upper RGB stars (IDenissenkov &: Herwiel |2004| ). It is not clear how thermohaline con- 
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vection can explain the phenomenon of Li-rich K giants given that its efficiency should be 
lower in the more rapidly rotating stars. These arguments have motivated our search for an 
alternative RGB mixing mechanism. 

In this paper, we use a simple model of toroidal magnetic field generation in a differen- 
tially rotating radiative zone of a bump-luminosity RGB star to obtain order-of-magnitude 
estimates demonstrating that the buoyant rise of magnetic flux rings being formed close to 
the local minimum in /i may be a good alternati ve to the ^He-driven thermohaline convection. 



A similar model has recently been proposed by iBusso et al.l (120071 ) (hereafter, referred to as 
BWNC). However, they assumed that a rising ring always stays in thermal equilibrium with 
its surrounding medium, and we argue that this leads to a substantial overestimate of the 
ring's radial velocity. We account for the impact of /i gradients and discuss the origin of the 
magnetic rings. We show that, as a mechanism for the RGB extra mixing, magnetic buoy- 
ancy has some advantages over thermohaline convection and, therefore, it is worth further 
investigating by means of multidimensional MHD simulations. 



The RGB Stellar Model 



Our background RGB stellar model, in the radiative zone of which the formation and 
buoyant rise of magnetic flux rings are studied, represents a typical metal-poor upper RGB 
star. It has the initial mass M = 0.8 Mq, helium and heavy-element mass fractions Y = 0.24 
and Z = 0.0005, and the luminosity log L/Lq = 2.085 corresponding to the age of 13.65 
Gyrs. The RGB extra mixing is free to work in this star because the H burning shell has 
already erased the H-proffie discont inuity in it. The model has been computed using the 
stellar evolution code described by iDenissenkov et al.l (120061 ). Its chemical element mass 
fraction, /i, and proffies immediately above the H shell are plotted in Fig. [TJ 

For extra mixing to dredge up material deficient in C but not enric hed in Na, as require d 
by the observed abundance patterns in the metal-poor field RGB stars (iGratton et al.ll2000l ). 
it has to reach a depth between 0.045 i?0 and 0.055 i?© (Fig. [T^). We will assume the 
mixin g depth = 0.05 -Rp (sh o wn by the vertical solid line in Fi g. [H). In previous works 
(e.g., Denissenkov fc Weiss 1996 : Denissenkov fc VandenBerg 20031 ). the mixing depth was 
specified using either the relative mass coordinate SM^^^ = (Mmix — ^c)/(^bce — ^c) or the 
logarithmic temperature difference AlogT = logT(rc) — logT(rjnix)! where the subscripts 
"c" and "bee" refer to the He core boundary and to the bottom of convective envelope, 
respectively. However, we have noticed that, when plotted as functions of r, the abundance 
profiles remain nearly stationary in spite of the slow mass infiow from the radiative zone 
that feeds the H burning shell. Therefore, we have decided to simply use the radius for the 
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RGB mixing depth specification. The value of r^i^ = 0.05 i?© corresponds to SM^i^ = 0.135 
and AlogT = 0.245. 

For subsequent estimates of various quantities characterizing the efficiencies of magnetic 
flux ring formation and buoyant rise in the radiative zone of our RGB model we need to 
know some of its structure parameters at r = rmix. These are summarized in Table [1] along 
with the parameter values at r = rbce- We have also listed the values used by BWNC. Note 
that they have considered a half-solar metallicity bump luminosity model with the initial 
mass of 1.5 Mq. 



3. Rotation in the Radiative Zone 



In our magnetic buoyancy model, it is assumed that a toroidal magnetic field in the 
radiative zone of an RGB star is generated by its differential rotation that stretches a pre- 
existing poloidal field around the rotation axis. To elaborate on the model, we therefore need 
an estimate of the degree of dif ferential rotation in t he radiative zone. For this, we use the 
rotation profile M2 presented by lPalacios et al.l (120061 ) in their Fig. 2 (the dotted curve in up- 
per panel C). It shows the angular velocity f2 as a function of 6M = {Mr — Mc) / (Mbce — Mc) 
in the radiative zone of a bump luminosity model with initial values of M, Y, and Z almost 
identical to ours. The model's rotational evolution was computed using the stellar evolu- 
tion code STAREVOL designed for ID simu lations of the angular mornentum and chemical 



eleme nt transport in stellar radiative zones (jSiess. Dufour. fc Forestinill2000l : iPalacios et al. 



20031 ) . The code takes proper account of the rotation- induced meridional circulation, hori- 
zontal and vertical turbulent diffusion, as well as atomic diffusion and mass-loss. The initial 
zero-age MS model for the M2 run has been assumed to be a slow solid-body rotator with 
the surface rotational velocity of 5 kms~^. This choice is a reasonable replacement for de- 
tailed computations of the pre-MS and early MS evolution during which a low-mass star 
experiences a strong magnetic breaking of its much faster initial surface rotation and a core- 
envelope rotational coupling leading to its nearly uniform internal rotation on a much shorter 
time scale than its MS life time. 



Another important assumption made by iPalacios et al.l (120061 ) is uniform specific angu- 
lar momentum distribution in the RGB convective en velope. This is required to explain the 
origin of rapidly rotating red horizontal branch stars (jSills fc Pinsonneaultll2000l ). This em- 
pirical demand for a strong differential rotation in the convective envelopes of low-mass RGB 
stars is supported by recent 3D hydr odynamic simulat i ons o f the interaction of turbulent 
convection and rotation performed by lPalacios fc BrunI (120061 ). 
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The M2 bump luminosity model has ^(rbce) ~ 10~^ rads~^, and fi(rinix) ~ 10~^ 
rads~^. We will use these values in our follo wing; discussion . So, in the absence of trans- 
port processes other than those considered by lPalacios et al.l ( 120061 ) . the angular velocity in 
the radiative zone of our RGB stellar model could vary with the radius approximately as 
fl{r) = fl{r^ix) {r^iy^/ r)^ . This steep rotation profile results from the conservation of angular 
momentum in the mass inflow and modest angular momentum redistribution by meridional 
circulation and vertical turbulent diffusion. 



4. Suppression of Thermohaline Convection by Horizontal Turbulence 



Denissenkov fc Pinsonneaultl (j2008bl ) have derived a general criterion for convective in- 
stability in the presence of a negative and strong horizontal turbulence, the latter being 
characterized by a diffusion coefficient ^ t'moi- They have allowed for a possibihty that a 
convective element of the diameter d can travel a distance / > d. One of the consequences of 
this criterion is that thermohaline convection may be suppressed by the horizontal turbulent 
diffusion unless 



Vrad - V 



1 

> -- 



ad 



3K + Dy 



(1^ 



wher e Vrad is the r adiative temperature gradient (compare this inequality with condition 5 
from I Vauclairl 120041 ) . Note th at a similar re sult can be obtained from the dispersion relation 
(10) derived and analysed by lUlrichI (119721 ). One should only take into account that, in the 
presence of the strong horizontal turbulence, Ulrich's molecular diffusion coefficient D has 
to be replaced with Dh- Then, for the thermally limited modes, one readily finds 



D 



the 



-Duirich X 1 



Dh Vrad — V, 



ad 



K 



IV,, 



(2) 



where -Duirich is Ulrich's original thermohaline diffusion coefficient. From the last expression it 
follows that thermohaline convection may operate {Dthc > 0) only when | V^|/(Vad — Vrad) > 
D^/K, which is close to the requirement ([T]), provided that < K. The difference in the 
right-hand-side ratios betwee r i the t wo conditions comes about from the fact that, unlike 
Denissenkov fc Pinsonneaultl (j2008bl ). Ulrich neglected the contribution of D to the heat 
diffusion. That was warranted for the nonrotating case considered by him, in which Z^h = 0, 
and therefore D = u^^i ^ K. 



Palacios et al.l (120061 ) have shown that Z^h always stays comparable to at r ^ r 
in their M2 RGB model. Moreover, the ratio D^^jK turns out to be larger than ~ 10 
at all radii between r^^^ and rtce, both in the bump luminosity model and in the more 



mix 
3 
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evolved model shown in their Fig. 5 (panels C and D). Given that the ^He burning can make 
the ratio |V^|/(Vrad ~ Vad) ~ 10~'^ at most, while the RGB extra mixing will necessarily 
reduce it below this limit (by smoothing out the /i-gradient), the linear analysis predicts 
that thermohaline convection may be suppressed in these models, especially in the vicinity 
of r-mix- It would be inconsistent to ignore this prediction because similar eroding effects of 
the horizontal turbulence on the model r otational and rnixing properties have already been 
included in the simulations performed by lPalacios et al.l (120061 ). They are responsible for the 
signif icant reduction of the efficiency of mixing by meridional circulation (jChaboyer fc Zahn 



19921 ) and for the erasing of the latit udinal differ ential rotation, the latter effect allowing to 



consider as a function of r alone (jZahij|l992l ). Besides, Zahn's concept of the rotation- 
induced anisotropic turbulence in stellar radiative zones, with horizontal components of the 
turbulent viscosity strongly dominating over those in the vertical direction, was repeatedly 
used in stellar evolution computations to facilitate the penetration of /i-gradi ent barriers by 
the vertical turbulent diffusion both in MS stars and in upper RGB stars (ITalon fc Zahn 



1997l : lTalon et all 119971 : lMaedei]l2003l : IPalacios et al.ll200l l2006h . 



5. Generation of Toroidal Magnetic Field 



Following iMestel &: Weisd (Il987l ). we assume that, like in the solar-type MS stars, the 
differential rotation in the radiative zone of the RGB star creates a toroidal magnetic field 
by shearing a pre-existing constant poloidal field Bp = {S^, Bq, 0} (we use the spherical 
polar coordinates). As a result of it s buoyancy-related undular instability (e.g., lAcheson 
19781 : ISpruit fc van Ballegooijeru Il982l ) , the toroidal field is prompted to form magnetic flux 
rings that will rise toward the bottom of convective envelope, thus producing chemical mixing 
and also participating in angular momentum redistribution. 



At the stellar equator {9 = 90°), the momentum and induction equations (jCharbonneau &: MacGregor 



19931 ) can be reduced to 



'dt 
db 

dt 



db 

dr ' 



(3) 
(4) 



where ujA,r = Br/ (y/Airpr) is the local Alfven frequency associated with the radial field 
component Br, and b = rB^/Br- In equations ([3]-|l]), we have omitted the viscosity and 
magnetic diffusivity, for these will be shown to work on much longer time scales than the 
formation and buoyant rise of magnetic rings. For the sake of simplicity, we have additionally 
assumed that the poloidal field's configuration is such that Be = at the equator. Neglecting 
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changes with the radius of the coefficients in equations ([3]-|l]), we find that the variations of 
Q and b can locally be described by the same Alfven wave equation 

So, at a given radius, the initial (continuing while t <^ ^^a^) decrease of and increase of \b\ 
can be approximated as fl{t,r) ^ Q^i^^^^r) cos{uA,rt), and b(t,r) ^ femaxl'^) sin(ciJA,r From 
([3]-|l]) we can also estimate the ratio of the wave amplitudes ^^max/^max ~ ^A,r/r, hence 
(-B<^)max ~ Br{flnia.x/^A,r)- This cusurcs that the sum of the rotational kinetic and toroidal 
field potential energy (per unit gram) is conserved, 

I (r^? + ^ = l (^^max)' . (5) 

Substituting the values of r = rmix and p = p (rmix) from our RGB model (Tabled]) into 
the above relations, we find that 



00^1 = 1.50x10^57^ yrrad"\ (6) 



(5^)max = 4.75X 10^(n^ax)-5 5r, (7) 

where (fimax)-5 = ^max/(10~^rads"^). Equation ([7]) gives an order-of-magnitude estimate 
of the toroidal field maximum strength that may be generated by the differential rotation 
at r = Tinix, while equation ([6]) estimates its growth time. For example, we may expect 
(from equation Hj) that, after a time At ^ ^aI' ^ toroidal field ^ Vt^^^B^q^t ^ 
(-B,^)max'^A,r At = 3.16 X 10^ (fimax) -5-Br 5' At wiU be Created, where q = {din fl/d\nr) < 
is the initial rotational shear, and At is expressed in years. Assuming that ilmax = ^('"mix) = 
10~^rads~^ (Tabled]) at the moment when the differential rotation begins to stretch the 
poloidal field, we obtain the estimate 

B^{At, r^ix) ^ 3.16 X lO^S^ q At, (8) 

where At <^ 1.50 x 10^ i?^^ yrs. The relative decrease of Q for the same period of time is 
|Afi/fii^ax| ~ 4.42 X 10-^52g2 (At)2 ^ 1, In the gene ral case of B f, ^ and ^] = n{r,9) 
we would have B^, ^ r sin 9 {Vn, Bp) At for At < ujI\ JSpruitlll999h . i.e. the both poloidal 



field components would be involved into the winding up of toroidal field. 

If the uniform rotation of the solar radiative core and the core-envelope rotational cou- 
pling in the low-mass MS stars are both produced by the back reaction of the azimuthal com- 
ponent of the Lorentz force emerging in the process of generati on of toroidal magnetic fields 



by the shearing of pre-existing poloidal fields, as proposed by ICharbonneau fc MacGregor 
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( Il993[ ). then the MS progenitors of low-mass RGB stars are required to possess poloidal 
magnetic fields with strengths of the order of 0.01 G to 10 G in their radiative interiors. The 
low magnetic diffusivity rj ~ 10^ sm^ s~^ in the low-mass MS stars excludes the ohmic dissi- 
pation of these fields. Therefore, the post-MS contraction of the H-exhausted core from its 
MS radius of ~O.2i?0 down to its RGB radius Tc ~ 0.02 i?© and magnetic fiux conservation 
might lead to the amplification of Br up to the values 1 G to IkG. From the same consid- 
erations, it also follows that B,. in the radiative zone of the RGB star may decrease with 
the radius as Br ~ -Br(^mix) ('"mix/"'")^- It is true that the mass infiow in the radiative zone 
may sweep the frozen-in poloidal field toward the H burning shell on a long time scale of the 
order of Ar/\r\ = {r^ce — ''^mix)/|'^| ~ 3 x 10^ yrs (Tabled]). However, the poloidal field may 
be replenished by a dynamo operating at the bottom of convective envelope. It may then be 
entrained and redistributed all over the radiative zone by the mass infiow, which could also 
produce the dependence Br oc r~^. These estimates will be used in our further analysis. 



6. The Buoyant Rise of Magnetic Flux Rings 



6.1. General Results 



Because we assume a differential rotation in the radiative zone, the inertial frame of 
reference has been chosen. The magnetic fiux rings are assumed to be axisymmetric with 
respect to the rotation axis. In the spherical polar coordinates their radial and latitudinal 
accelerations are 



dUr 

~dt 



dug 

~dt 



+ 



u 



r 



sin^ e + 



Pc 



GMr 



Pc + P \ 



n^r sin^ 9 



/ 9 I 2 
Ur\/Uf. + Ui, 



47rr(pc + p) Tca pc + p 



+ — cot 6 — \ rr sm 6 cos 6 



(9) 



pc- P 
Pc + P 



n'^r sin 6 cos 9 



Airripc + p) 



cot I 



Cp Pc 
ira pc + p 



U0\ + ui. 



(10) 



Thes e equations (jChoudhuri fc Gilmarull987l : lMacGregor fc Cassinellill2003l : lMacDonald fc MuUan 
20041 ) take into account the centrifugal reduction of the local gravitational acceleration, the 
buoyant force, the magnetic tension force, and the aerodynamic drag force (we employ the 
drag coefficient Cd = 1). They are supplemented with the equation 



^{u^r sin 6*) = 
dV ^ ' 
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describing the conservation of the azimuthal component of the specific angular momentum. 
The subscript "e" means that the respective quantity is referred to the external medium 
surrounding the rings, while all other quantities, except Q, refer to ring properties. 

At its starting position (ro, 6^0), a ring is specified by its initial cross-section radius 
ao, the strength of the frozen-in toroidal magnetic field -B^,0; and internal thermodynamic 
properties: Pq, Tq, po, and /xq- The dynamic equilibrium between the ring and its surrounding 
medium requires that 



Bl AP (P^-P) 1 1 



where j3 = P/(-B^/87r) ^ 1 is the ratio of the thermodynamic to magnetic pressure in the 
ring. The conservation of mass and magnetic flux of the ring determine how its radius and 
toroidal field evolve during its motion 



a\ _ Poro sin 6*0 _ fi^,o 



ao / pr sin 9 B, 



(13) 



where vq < r < rbcc- For the equation of state, we use the ideal gas law. Finally, the entropy 
change in the ring is described by the simple equation 

f^?^^l^^ST^2,K% (14) 



derived bv lMacGregor fc Cassinellil (120031 ). Here, S = In [(P/Po)(Po/p)^], ST = {T^-T)/T^, 
and 7 = 5/3. 

The initial conditions for the ring are specified by ro = r^iy^ = 0.05 Rq, < Oq < 90°, 
Ur = 0,ug = 0, = rosm6ofl{ro, 60), (3 = Po, a = do, S = 0, 6P = {Pe-P)/Pe = (l+/5o)"S 
6p = {pe—p)/pe = {^+Po)~^, and ST = 0. The ring's radius a will be measured in units of the 
local pressure scale height Hp = P^^/gp^, where g = GM^/r'^ is the gravitational acceleration. 
The implicit assumption of uniform pressure inside the ring requires that a <^ Hp (the thin 
ring approximation). For brevity, we will denote a = {a /Hp) and, given the previous remark, 
we will only consider cases with a <C 1. 

The above differential and algebraic equations have been solved numerically. The motion 
of the ring in the equatorial plane (^0 = 90°) is shown in Fig. [2] for our assumed rotation 
law Q{r,9) = ^{r^ix){r^ix/ry, in which case 9(t) = 9q. Solid curves correspond to the 
initial value of ao = 10"'^, while dashed curves have ao = 10^^. The powers of ten near the 
curves give their specified values of Pq, the respective toroidal field strengths (in MG) being 
displayed in parentheses. Our computations have shown that for /^o ^ 10^ (in other words. 
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for (-B^)o ^0.1 MG) and ao > 2 x 10 ^ the dependence of the total buoyant rising time on 
Po and ao can be approximated as 

- 1.4 X 10^ f (J^y y, (15) 



^108 J VlO 

This is transformed into the average buoyant velocity 

{vh) = — = ^ ~ 15 — cms , (16) 



which is nearly five orders of magnitude smaller than the appropriately scaled average velocity 
used by BWNC (for the scaling, we have used the data from the fifth column of Table [1] and 
the column RGB-1 of Table Ej). We have tested that this big difference is entirely caused by 
the assumption of thermal equilibrium between the ring and its surroundings made in the 
cited paper. 

The dashed and dot-dashed curves in Fig. [3] show the dependences of tb on ao for 
jSo = 10^ and /?o = 10'' obtained by solving equations (P- fTTi) . while the dotted lines represent 
their approximations by ( ITSi) . It is seen that the exact solutions strongly deviate from 
the approximate ones at Oq ^ 2 x lO"'^. This limit corresponds to a regime in which the 
aerodynamic drag force comes into play. Indeed, the last terms on the right-hand-sides of 
equations (P-[TOl). that describe the drag force, are negative and inversely proportional to a. 
They lead to a slowing down of the ring's buoyant rise at small a, which is not reflected in 
equation (ITSi) . 



6.2. Rings with ji Reduced by ^He Burning 

Our numerical computations can optionally take into account the fact that the rings 
carrying the nuclear processed material are necessarily formed in the region of the local /x 
depression maintained by the ^He burning (Fig. Wp)- They must therefore have a lower n 
than the bulk of the radiative zone through which they move. As a result, their average 
buoyant velocity is found to weakly depend on the toroidal field strength. However, the 
frozen-in toroidal field is still needed for them to remain cohesive while rising. 

Dotted curve in Fig. [2] shows the path of a ring with Afi = fi^ — fJ' = 5 x 10~^, = 10""^, 
and Po = 10^° that corresponds to (-B^)o ~ H kG. Its buoyant rising time is ~ 1.3 yrs. 
Approximately the same short time is obtained for a ring with (-B^)o ~ 1-1 kG {f3o = 10^^). If 
those two rings had = /le then their rising times would be much longer and vastly different, 
namely: tb ~ 10^ yrs, and tb ~ 10^ yrs, respectively (eq. [T^ . This means that their buoyancy 
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is now controlled by the jj, difference alone. Let us neglect, for a moment, the difference in 
temperature between the surrounding medium and the ring, 6T = (Tp — T)/Tc ~ 0. The 
buoyant acceleration is Ob ~ g5p. In the case of A/x = 0, the initial positive difference in 
density 5p ~ 5P ~ (3^^ is simply due to the excess magnetic pressure inside the ring (eq. 
[T2!) . On the other hand, in the case of A/x > 0, and for a weak toroidal field {5P ~ 0), 
we have 6p ^ 6p = Ap/pe- We can introduce the effective toroidal magnetic field B^s 
associated with a specified value of Sp such that Sp = P~g. It turns out that -Befr ~ 9.8 
MG for Ap = 5 X 10"^ at r = rmix in our RGB model. Hence, toroidal magnetic fields with 
{Bip)o ^ Bcs will have a negligible effect on the motion of the ring with the reduced p. For 
Ap = 5 X 10~^, /5o = 10^°, and Qq = 10^^ the rising time tb ~ 1-6 yrs is still short. However, 
it increases up to tb ~ 33 yrs and tb ~ 132 yrs for rings with the initial radii oq = 5 x 10~^ 
and ao = 10~^, respectively, following the dependence tb oc Oq for the thicker rings (eq. [T51) . 

The role that the reduced p plays in the acceleration of the buoyant rise of the ring 
containing a weak magnetic field can be elucidated if we consider, for simplicity, that the 
ring's vertical motion consists of the following recurring sequence. Initially, let the ring have 
6T = but 6p > 0, hence 6p = 6p + Pq-^ > 0. Now, let the ring rise adiabatically until 
the accumulated difference in T compensates that in p, i.e. ST = 6p. At this moment, the 
ring stops and waits a while for the heat exchange to make 6T = before starting to rise 
adiabatically again, and so on. But the waiting time is in fact the thermal time rth, therefore 
it is inversely proportional to the heating rate dS/dt that linearly depends on the ratio bT jc? 
(eq. [H]). When the ring stops, it has 5T = 5/i + /3q"^ Therefore, if bp ^ Pq-^ then 6T ^ 6p, 
hence we have tb oc Tth oc 6p^^ = /5efj a^. On the other hand, if 6p <^ /Sq-^ then 6T ^ P^-^, 
and tb oc rth oc Poo^ (eq. [15]). The transition from the one to the other regime occurs at 

p^' ^6p = p;l 

We have checked that equation (ITH) gives a correct order-of-magnitude estimate (this 
time, at >8x 10~^ though) for the buoyant rising time of a ring with a reduced p provided 
that Po is replaced by Pes = Sp~^ (the dotted line approximating the solid curve in Fig. [3]). 
It is also important to note that the coefficients in equations f|T5] - [T6l) have been obtained 
for a particular RGB model (our model from Tabled]), therefore they are model dependent. 
For example, given that tb is expected to be inversely proportional to the thermal diffusivity 
K as well (eq. [H]), the latter being roughly proportional to the luminosity, we predict that 
these coefficients should change by a factor of 10 (tb decreases, while (fb) increases) toward 
the RGB tip (log L/Lq ~ 3.3). This scaling is indeed confirmed by our computations. 
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6.3. Comparison with Results Obtained by BWNC 



Given that one rising magnetic flux ring carries the mass mb = 27r^ro po sin = 
27r^ra^psin ^, N such rings present in the radiative zone at the same time will provide 
chemical mixing with the mass rate 



NirtY 



(17) 



where tb is their buoyant rising time. The quantity Mb has to match the observationally 
constrained rate of the RGB extra mixing Mmix ~ 4 x lO^^Mgyr"^ (BWNC). The right- 
hand-side of equation ( ITTl) is a function of the ring parameters A^, oq, and {B^)q. For a 
ring with a reduced /i, the third parameter should be replaced with A/x, unless -C l//5o 
(§ 16. 2p . In the case of A/i = 0, considered by BWNC, the equating of Mb to Mmix gives a 



relationship between A^, oq, and (-B^)o- To estimate a reasonable value of ao at r 



they 



have referred to the characteristic dimension a© ~ 1000-2000 km of magnetic flux tubes that 
are believed to exist deep in the solar convective zone. Beside the RGB model, they have also 
considered a representative model for low-mass asymptotic giant branch (AGB) stars, whose 
structure above the H burning shell resembles that of upper RGB stars. There is indirect 
evidence, such as distinctive ^^O/^^O, ^^O/^ ^O, ^^C/^^C, and jM/C abundance ratios in the 



meteorite grains of AGB circumstellar origin (|Wasserburg et al. 



of carbon-enhanced metal-poor stars (IRvan et al 



2006: Denissenkov fc Pinsonneault 2008a: Lebzelter et al. 



2005 



20061) and in the atmospheres 



Sivarani et al.ll2006l : iMasseron et al. 
2008h. indicating that extra mixin g 



may operate in the radiative zones of these stars too (INollett. Busso. fc Wasserburg||2003l ). 
In their RGB-1 and RGB-2 magnetic buoyancy models, BWNC have employed the values 
of obtained assuming that a(rbce) = clq using equation (IT^il with the stellar structure 
parameters from their AGB and RGB stellar models, respectively. These values are listed in 
our Table [2] along with the corresponding estimates for our RGB model. 

BWNC have also assumed that there is only one magnetic flux ring floating in the 
radiative zone at any time. Substituting the number A^ = 1 together with the values of 
from Table [2] and M^ix = 4 x 10~^ Mq yr~^ into equation f[T7l) constrains the required 
values for ty, and (fb) (Table [2]). Under the assumptions, made by BWNC, that the ring 
always stays in the thermal equilibrium with its surrounding medium {6T = 0) and that 
the only force impeding its motion is the aerodynamic drag force, it is easy to show that 
tb oc a/Z^o /flo oc (5^)0/1/00. It is from the last relation that BWNC have determined the 
toroidal fleld strengths needed to drive the RGB extra mixing by magnetic buoyancy (Table [21 
{B^)o is transformed to (i?i^)bce using equation [13]). They pointed out that these results are 
consi stent with existing observations of magnetic flelds in red giants (e.g., 
200lh . 



Blackman et al. 



- 15 - 



However, as we have noted, BWNC underestimated the ring's rising time by assuming 
that the heat exchange between the surrounding medium and the ring occurs instantaneously, 
which maintains AT = all the time. Our more conservative assumption explicitly takes into 
consideration the radiative heat exchange, which leads to a much longer ring's rising time, 
as approximated by equation (fTSl) . Employing the values of oq, tb, and Pq from the columns 
RGB-1 and RGB-2 of Table [21 we find that the rising times have been underestimated by the 
factors 6 x 10^ and 4 x 10^ for these buoyancy models. To keep the same estimates for the 
toroidal field strength, the ring number N has to be increased by the corresponding factors. 
However, the volume occupied by ~ 10^-10^ rings is comparable to the total volume of the 
radiative zone, in which case the stellar structure would be greatly disturbed, especially at 
the bottom of convective envelope where the rings encounter turbulent convection. Here, 
the frozen-in toroidal field is either quickly dissipated via the strong turbulent diffusion, thus 
depositing its energy at the radiative/convective interface, or it inhibits c onvective motions 
if its potential energy exceeds the turbulent kinetic energy (IMossi l2003l ). In either case, 
the stellar structure would be strongly modified. In the next section, we will show that a 
reasonably small number of N can be obtained only for rings with a reduced /x. 



6.4. There Still May Be a Solution 



By assigning = 1, BWNC have implicitly assumed that the average time tfomi needed 
to form a magnetic flux ring is equal to its buoyant rising time ty,. Indeed, instead of (fTTll 
the buoyancy mixing mass rate ought to be calculated as Mb = mh/tf„rm = Nrrihlth, where 
N = tb/i^form- We assume that, hke in the case of the solar tachocline (e.g.. lSchmitt fc Rosner 



19831 ). the appropriate MHD mechanism responsible for the formation of magnetic flux rings 
in the vicinity of the H burning shell is the undular buoyancy instability. The crite r ion fo r 
its development has been extens i vely discussed in t he literatur e, e.g. by lAchesonl (119781 ) . 
Spruit fc van BallegooijenI (119821 ) . ISpruitI (Il999l ). and iFanI (120011 ). It has been shown that a 
diffusive toroidal magnetic field gets broken into distinct arching flux tubes when 



ip ) crit 



r K 



provided that {d\nB^/ d\n.r) = 0(1). A profile of (-B^)crit in the radiative zone of our RGB 
model is plotted with the dashed curve in Fig. UK. Its corresponding profile of jScrit is shown 
with the dot-dashed curve. At r = r^i^ (log r^i^/ Rq ~ —1.30), the critical toroidal field is 
254 kG and P^nt = 1-8 x 10^ (Table [2]). 

In the presence of rapid rotation, such that Q ^ iJA,<fi = B,p/( x/4:7ipr), the Coriolis 



force reduces the instability's growth rate ujA,ip by the factor ua^^/Q ( jPitts &: Taylerlll985 
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SpruitI Il999l ). In our RGB model, fi(rmix) = 10 ^rads ^ strongly exceeds cja. 



2.11 X 10~^^ -B,^ rads^^ for all reasonable strengths of ^ 47 MG. Therefore, we use the 
estimate of tinst ~ ^/^\ip for the instability growth time. The winding up of a toroidal 
field by differential rotation continues until -B,^ ~ QBrAtq reaches the critical value f[T5]) . 
This takes (At)crit ~ 3.17 x 10~^ (i?<^)crit/(-Br ^-5 ^) years. After that, it will take another 
^inst ~ 1-93 X 10^° [p {r / RqY Q_^] /{B^)^^.^^ years for the buoyancy instability to occur and 
assemble the rings. So, the total ring formation time is tform = ('^'^)crit + ^inst- For our 
simple estimates to be true, tform should be shorter than u!^\. = 7.82 x 10'^ ^/p {r / Rq) / Br 
yrs, which allows to consider that the rotational shear is nearly constant and that, for our 
assumed rotation law in the radiative zone fi_5 = 10^ ('"mix/'")^; Q = 0(1). For this case, and 
assuming that B^ = 7.32 (r^i^/r)^ G, the four characteristic time scales are plotted in Fig.Hjo. 
In particular, we have tform = 2.20 yrs at r = r^i^. Note that (At)crit oc {Brfl)~^, tinst ^, 
and u!^\ oc B~^. We have used the parameter Br{rmix) = 7.32 G for which the time scales 
(At)crit and tjnst coincide. Its value scales as Q~^{r^i-^). At a fixed value of Q{r^i^), -Br('^mix) 
determines which of the two time scales, (A)f:)crit or tinst, niakes a predominant contribution 

to tform- 

It is very likely that the buoyancy instability forms not just one but a number n > 1 
of magnetic flux rings during the time tform, therefore = ^(tb/iform)- The equating of 
Mb = nrrib/tfoTra to Mmix = 4 X 10~^ Mq yr~^ constrains the ring's minimum radius as a 
function of tform, n, and the initial colatitude 

' ^ M , (—Y = 9.78 X 10"^ (19) 

Hp J n 27r^rQpo sm t^o V '"o / n suiOq 

where tq = rmix, and tform is expressed in years. 

For the values of tform = 2.20 yrs, n = 1, and 9q = 90°, we calculate Oq = 1.47 x 10~^ (in 
units of Hp, as usually), and the total number of rings in the radiative zone A^ = 0.454 tb- As 
we have mentioned before, the ring with the radius oq = 10~^ and p reduced by Ap = 5x 10^^ 
has the buoyant rising time tb ~ 132 yrs. This means that the ring with ao = 1.47 x 10^^ and 
the same mean molecular weight would cross the radiative zone in ~ (1.47)^ x 132 = 285 
yrs (eq. [T^ . Hence, we have A^ = 129. Note that this number does not change if the values 
of the parameters tfomi and n are varied, provided that our buoyancy mixing mass rate is still 
constrained to match Mmix. Indeed, in this case both the ratio tform/n (eq. [19]) and (eq. 
[TSll are proportional to which leads to its cancellation in the relationship A^ = r;,(tb/iform)- 



We have shown that the reduced mean molecular weight in the magnetic flux rings 
formed in the region of '^He burning accelerates their buoyant rise quickly enough for their 
total number needed to maintain the RGB extra mixing to be reasonably small (A^ ~ 10^ 
rings would only occupy ~ 10^^-10"^ part of the radiative zone). However, the estimated 
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value of ~ 1.3 x 10^ is in fact a lower limit obtained under the most favourable assumptions. 
In particular, we have silently assumed that all rings have the same formation time and cross- 
section radius related by equation (HM . and that no rings are formed above the radius r^i^. 
These assumptions are obviously not true. Relaxing either of them will certainly increase 
the total ring number in the radiative zone. To figure out what effect may be produced by 
relaxing the first assumption we would have to determine a relationship between tform and 
Gq for a spectrum of ri ngs created by the buoya ncy instability, e.g. like it has been done for 



the solar tachochne by lSchmitt fc Rosneii (119831 ). This problem is out of scope of the present 
preliminary study. 

We find it more important to address here the second issue. Indeed, given that the 
critical strength of toroidal magnetic field for triggering the buoyancy instability decreases 
rapidly with the radius (dashed curve in Fig. Hti), while our estimated ring formation time 
stays much shorter than the toroidal field growth time uj^\ (Fig- Hb), the formation of 
magnetic flux rings at vq > r^aix appears to be unavoidable. These rings would contribute 
to the total ring number present in the radiative zone at the same time but they would 
not participate in the chemical element transport because their constituent material has not 
been nuclearly processed. In fact, the buoyancy of these "parasitic" rings should be reduced 
compared to the buoyancy of the rings originating at r ^ rmix because their mean molecular 
weight does not differ from that of the surrounding medium through which they rise. On 
the other hand, they are formed in a region where the thermal diffusivity K is higher than 
at r = rjnix, therefore their heat exchange with the surroundings goes faster, which should 
accelerate their buoyant rise. Our computations show that the rings formed in the region 
of constant /i (at r > rmix + 0.05 i?0) rise 10^ to 10^ times slower than the rings formed in 
the /x-depression domain. Given that their formation takes between 10 and 100 yrs (solid 
curve in Fig. |Dd), a number of the "parasitic" rings present in the radiative zone at the same 
time may be as large as 10^. Unfortunately, our model is too simple to be able to predict if 
and how the "parasitic" rings could impede the large-scale magnetic buoyancy mixing, but 
we do realize that they may create a real problem for the mechanism proposed by BWNC 
that has been revised in our paper. Indeed, although BWNC did not discuss a formation of 
magnetic rings, it is difficult to understand why the rings can not be formed at r > r^^^ as 
easily as they are created at r = rmix- 



7. Discussion 



Our approximate analysis of the magnetic flux ring formation and buoyant rise in the 
differentially rotating radiative zone of the bump-luminosity RGB star gives some support 
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to a combined "magneto-thermohaline" mode of the RGB extra mixin g, as opposed to the 
pure th ermohahne and pure magnetic buoyancy modes proposed by ICharbonnel fc Zahn 
( l2007al ) and by BWNC, respectively. A key component to the operation of our mixing 
mechanism, which is certainly present in all upper RGB and low-mass AGB stars, is the 
mass inflow in the radiative zone. It makes two important things. First, the conservation 
of the specific angular momentum in the mass inflow results in a steep f2-profile with a 
rotational shear q ^ 0(1) that do es not appear to be strongly reduced by the m eridional 
circulation and turbulent diffusion (jPenissenkov fc Toutll2000l : iPalacios et al.ll2006l ). Second, 
the mass inflow entrains a diffusive magnetic field that is pro bably generated via a.n a-Q or 
dynamo just beneath the bottom of convecti ve envelope (INordhaus. Blackman. &: Frank 
20071 : iNordhaus fc Blackmanll2007l : lMestelll200ll ). A very high ratio of the ohmic dissipation 
time to the average inflow time ((Ar)^/ {ri))/{Ar/ (|r|)) ~ 1.2 x 10^ guarantees that this field 
will be redistributed over the whole radiative zone. Because the volume occupied by this field 
is squeezed when the flow approaches the H burning shell, the field may become stronger 
at a smaller radius. In general, thi s random fie ld has an unstable configuration that will 
decay on a short Alfven t ime s cale (jSpruitI 119991 ). However, in their 3D MHD simulations 
Braithwaite fc Nordlundl (120061 ) have shown that such unstable random field may evolve 
into a stable "twisted torus" configuration with toroidal and poloidal field components of 
comparable strength. The poloidal field lines get wrapped around an axisymmetric torus, 
thus forming an approximate dipole. These simulations have provided the first plausible 
explanation of the origin and stability of dipole magnetic fields in presently non-convective 
stars. It is assumed that the seed random magnetic field in those stars had been a dynamo 
left over from the period of their protostellar convective contraction. 

Our RGB extra mixing model has much in common with the solar magnetic spin- 
down mod^]_that_20^sti^^ field in the solar radiative 
core (IMestel fc WeissI 119871 : ICharbonneau fc MacGregorl Il993l ). Note that even the values 
of T-bce = 0.996 i?0 and ^(rbce) = 10~^rads~^ used in our RGB model are close to the cor- 
responding solar values. The important difference is the assumption of strong differential 
rotation in the radiative zone of our model, while th e present-day Sun is known to be a 
nearly solid-body rotator, at least above r = 0.2 Rq JCouvidat et al.lboosh . However, the 
young Sun did possess a strong differential rotation in the core which had resulted from its 
spinning up during the pre-MS contraction and angular momentum loss from the surface 
via a magnetized stellar wind. The magnetic spindown model assumes that the differen- 
tial rotation in the young Sun was broken by the back reaction of the Lorentz force that 
emerged when the differential rotation was windi ng up a toroidal field from the pre-existing 
poloidal field (ICharbonneau fc MacGregon 119931 ). or by the magneto-rotational instability 



jMenou fc Le Mei]l2006h 
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Whereas an important role in damping large-s cale toroidal field oscillations in th e young 
Sun is thought to be played by the phase mixing (ICharbonneau &: MacGrego ] |l993h . we do 



not think that this is also true for our magnetic buoyancy model. Indeed, the phase mixing 
time scale at r = rmix in our RGB stellar model is 

*p = 2 2 ~ 2.9 X 10^ yrs, (20) 



where we have used an estimate of tp obtained by ISpruitI (Il999l ). assuming that qa = I- This 
is comparable with the time needed for the whole radiative zone to be thoroughly mixed, 
^mix = (Afbce — Afmix)/A!fmix ~ 5.5 X 10^ yrs. Heucc, long before the toroidal field oscillations 
on neighboring magnetic surfaces get out of phase the magnetic buoyancy instability will 
come into play, the magnetic fiux rings will be formed and reach the convective envelope. 

Compared with the pure thermohaline mixing, our model has the following advantages. 
First, the horizontal turbulent diffusion is unlikely to hinder the buoyant rise of magnetic 
flux rings with a reduced n because the frozen-in toroidal fleld will not allow turbulence to 
penetrate plasma in the rings and decrease the /i contrast between the ring and surrounding 
material. On the contrary, there is nothing to prevent the horizontal turbulent diffusion 
from eroding the /i contrast in thermohaline convective elements (§|4]). 

Second, thermohaline convection may fail to explain the operation of enhanced extra 
mixing in rapidly rotating Li-rich K-giants, which is needed to activate the Cameron-Fowler 
mechanism, be cause one would expect that the growth of "salt flngers" is impeded by rotation 



(ICanutd Il999l ). Oppositely, it would be natural to suppose that magnetic flux rings are 
formed more efficiently in the more rapidly rotating stars. This hypothesis is supported by 
the fact that the fastest rotators among young cluster solar-type stars appear to have the 



short est time scale of rotational coupling between the core and envelope (e.g., Ilrwin et al. 



20071). 



Third, whereas the efficiency of thermohaline convection in upper RGB stars is depen- 
dent exclusively on the abundance of ^He left in the radiative zone, magnetic buoyancy can 
in principle be driven by differential rotation alone, provided that it succeeds in winding 
up a sufficiently strong toroidal field. By the end of the RGB evolutio n of a low-mass star. 



its eny elope He abundance gets depleted by a factor of ten or more fICharbonnel fc Zahn 



2007al ). So, what will then drive extra mixing in this star on the AGB? The higher rate 
of heat exchange between "salt fingers" and their surrounding medium increased in propor- 
tion to the luminosity will not help because the evolutionary time scale decreases inversely 
proportional to the luminosity which cancels the former effect. Thus, there is a need for 
an additional driving parameter that would not let extra mixing die out. Such parameter 
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might be the rotational velocity. Indeed, it is known that a surprisingly fast rot ation has 



somehow survived in red horizontal branch stars , in spite of the RGB mass loss ( iPeterson 



19831 : ISills fc Pinsonneaultll2000l : iBehr et al.ll2003l ). It could be used on the subsequent AGE 
evolutionary phase to drive our magneto-thermohaline mixing. The much lower (compared 
to the RGB phase) ^He abundance left in these stars could be compensated by a stronger 
toroidal magnetic field or a larger number of magnetic flux rings generated in the presence 
of a higher mass inflow rate (Mmix ~ lO~^M0yr~^, according to BWNC). 

Fourth, thermohaline convection cannot penetrate below the radius rmin ~ 0.063 i?0 at 
which n has a minimum (the right vertical dotted line in Fig. [1]). However, mixing down 
to this depth would be t oo shallow to reduce the surface carbon ab undance (Fig. Wp), as 
required by observations ( Gratton et al. 2000; Smith &: Martell 20031 ). An overshooting on 
a length scale of order Hp could solve the problem (the left vertical dotted line in Fig. [1] is 
placed at a distance Hp below r^in) but then thermohaline "fingers" would have to penetrate 
a region of higher /i where they experience a strong breaking. It should also be noted that the 
penetration of a region with the positive below r^i^ would reduce the average mixing rate 
-Dthc oc |V^| by decreasing the slope of the negative in the mixed radiative zone. In the 
magnetic buoyancy model, there is at least a potential possibility to dredge up the nuclear 
processed material from below r^[^. This can be done by rings with a frozen-in magnetic 
field of a few MG if such are formed closer to the major H burning shell. Of course, all the 
above assumptions need to be verified by more rigorous models. 

Fifth, it is interesting to note that the operation of our RGB mixing mechanism seems 
to produce environme nt in which the funct i oning o f the ^He-driven thermohaline convection 
is impossible. Indeed, ICharbonnel fc ZahnI (j2007bj) have estimated that a toroidal magnetic 
field of ^ 100 kG would entirely suppress the thermohaline mixing at r '"mix, while a 
field of 10 kG would be sufficient to inhibit the thermohaline instability in the upper one 
third of the narrow /i-depression domain. These fields are weaker than those needed for our 
mechanism to work. 

It is also important to note that the fiux rings considered in our paper are a convenient 

proxy for the more realistic f2-shaped loops that are the likely prod ucts of the magnetic buoy- 

ancy instability that operates in a layer of strong toroidal field (e.g.. lCaligari. Moreno-Insertis. &: Schussler 

19951 ). Downfiows/upfiows that take place in the loop legs during formation and rise may 

affect the efficiency of transport of processed material relative to the behavior obtained in the 

case of a flux ring. It is not clear how the strong rotational shear that is present within the 

radiative zone would affect this process. In principle, the only place where significant shear 

is required is within the layers where n is depressed and fiux tube formation takes place. If 

the overlying portion of the radiative zone were in a state of near-uniform rotation, the loop 
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legs might form a kind of conduit connecting the region containing processed material with 
the bottom of the convection zone. In any case, it is difficult to give an air-tight argument 
for what the rotational state of the radiative zone should be, either differential rotation of 
the kind assumed in our paper or near-uniform rotation enforced by the fields, that permeate 
the region. 

Another aspect of our computations that should be noted is that the results of sections 
6.1 and 6.2 were obtained for rings initially in thermal equilbrium with their surroundings, 
a state of maximal initial buoyancy. If, alternatively, the rings were in an initial mechanical 
equilibrium state (zero net force), the rise would be slower by virtue of the heating required 
to overcome the neutral buoyancy at t = 0. This would also increase the total ring number. 



8. Conclusion 



In this work, we have presented a simple model of the formation and buoyant rise of mag- 



netic fiux rings in the radiative zone o 
on ide a s and equations publis he d by 



( 1983 ). iMestel fc WeissI fll987l). IChoudhuri fc GilmanI ^ 987). ICharbonneau fc MacGregoi 



the bump luminosity RGB star. Our model is based 



19931) ■ ISpruitI (ll999l ). lEggleton. Dearborn. &: Lattanzid (120061 ) . BWNC, and lDenissenkov &: Pinsonneault 



Spruit fc van BallegooiienI (119821). ISchmitt &: Rosnei 



( l2008bl ). It qualitatively describes a possible mechanism for the RGB extra mixing, which 
we call the magneto-thermohaline mixing, as an alter native to the pure ^He-driven thermo- 



haline convection that has recently been proposed by ICharbonnel fc ZahnI (j2007al ). For our 
mechanism to work, the radiative zone has to possess a strong differential rotation and a 
poloidal magnetic field Bp > 1-10 G. We assume that the differential rotation stretches 
the poloidal field around the rotation axis, thus creating a strong toroidal magnetic field 
^ 0.1-1 MG. When the latter exceeds a critical value, the buoyancy-related undular 
instability comes into play to form magnetic fiux rings. These rings turn out to be buoyant, 
therefore they rise toward the bottom of convective envelope. 

We have shown that, when the radiative heat exchange between the ring and its sur- 
rounding medium is taken into account, the ring's buoyant rising time increases by about 
five orders of magnitude compared to the case considered by BWNC, when the ring and 
its surrounding medium are assumed to be in thermal equilibrium all the time. However, 
given that our model still neglects possible internal heating of the ring's material by resid- 
ual nuclear reactions and anisotropic thermal exchanges in the presence of strong oriented 
magnetic field, while it uses the aerodynamic drag coefficient that is 20 times as large as the 
one employed by BWNC, it is fair to say that BWNC might have fixed a safe upper limit 
while our paper fixes a conservative lower limit for the magnetic ring's rising velocity. 
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We have found that the number of rings needed to be present in the radiative zone at 
the same time to produce the observationally constrained rate of the RGB extra mixing is 
unreahstically large unless these rings originate from the region of the fi inversion maintained 
by the ^He burning. Such rings have a deficit of the mean molecular weight compared to the 
bulk of the radiative zone through which they move. Their buoyancy is mainly caused by the 
difference in n rather than by a deficit in density due to the excess magnetic pressure. The 
frozen-in toroidal magnetic field is still needed for the rings to remain cohesive while rising. 
That is why we have coined the term "magneto-thermohaline" mixing. Our model has some 
advantages over the pure thermohaline mixing model, the most important of which being 
the robustness of the magnetic rings against the eroding effect produced by the horizontal 
turbulent diffusion. Leaving aside the problem of the "parasitic" rings that are formed at 
^ ^ ''"mix + 0.05 Rq, our model looks promising. However, because it is based on a number of 
assumptions whose legality is impossible to confirm in the framework of our ID computations 
we call for its future verification by 3D MHD simulations. 

PAD and MP acknowledge support from the NASA grant NNG05 GG20G. The National 
Center for Atmospheric Research is sponsored by the National Science Foundation. 
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Fig. 1. — (a) The element mass fractions, (b) the mean molecular weight, and (c) its loga- 
rithmic (with respect to the pressure) gradient as functions of the radius in the vicinity of 
the H burning shell in our bump-luminosity RGB model. The vertical solid line shows the 
observationally constrained depth of the RGB extra mixing rmix = 0.05 i?©, the right dotted 
line - the radius rmin of the minimum /i, while the left dotted line is placed at r = rmin — Hp. 
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Fig. 2. — The motion of the equatorial magnetic flux rings with the cross-section radii 
flo = 10^^ i/p (sohd and dotted curves) and = 10^'^ Hp (dashed curves). The powers of 
ten show the specifled values of /So, while the numbers in parentheses give the corresponding 
strengths (in MG) of the frozen-in toroidal magnetic field. The dotted curve corresponds to 
the ring with = 10^° {B^ ^11 kG) whose mean molecular weight has been reduced by 
A/i = 5 X 10-^ 
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Fig. 3. — The dashed and dot-dashed curves show the dependences of the buoyant rising 
time on the ring cross-section radius for the rings with i3q = 10^ and i3q = 10^, respectively, 
whose chemical composition has not been changed by nuclear reactions (A/i = /ie — /i = 0), 
while the solid curve corresponds to the ring with /3o = 10^'' and A/i = 5 x 10^^ (/i reduced 
by ^He burning). These three curves represent our numerical solutions of equations (P-fTTl). 
The dotted lines give their corresponding approximations by equation (fTSl) . in which has 
been replaced by PeS = = (/i/A/i) for the third case. 
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Fig. 4. — (a) The critical toroidal field (dashed curve) and parameter (3 (dot-dashed curve) 
for the triggering of the buoyancy instability (equation [18]) , and the evolution of the toroidal 
field strength in the rising ring (solid curve, equation [T3|) . (b) The time scales for the winding 
up of the critical toroidal filed (Atcrit) and for the development of the buoyancy instability 
(^inst), as well as the total ring formation time (tform = Atcrit + ^inst)- These time scales 
should be shorter than the toroidal field growth time ojj^^.. In these computations, it has 
been assumed that = 7.32 (fmix/r)^ G, which makes (At)crit = ^inst at r 
log(rmix/^0) ~ -1.30. 



r^ix, where 
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Table 1. RGB Model Structure Parameters 


Parameter 


Units 




Our Model 


BWNC Model 


rmix = 0.0500 Rq rbce = 0.996 Rq 


rmix = 0.0495 Rq rbce = 0.912 Rq 


Mr 


Mq 


0.309 


0.331 


0.252 0.256 


Hp 


r 


0.257 


0.449 


0.371 0.622 




AMU 


0.593406 


0.593474 




T 


K 


2.16 X 10^ 


2.02 X 10^ 


2.51 X 10^ 2.26 X 10^ 


P 




14.8 


1.03 X 10-2 


5.19 2.55 X 10-3 


P 


dyncm-^ 


4.49 X 10^^ 


2.93 X 10^2 


1.87 X lO^*^ 8.49 X 10" 


K 


cm^ s 


9.35 X 10^ 


8.15 X 10^° 




ATI 




5.40 X 10"^ 


1.13 X lO-'' 








1.48 X 10"^ 


0.00 






2 —1 

cm s 


1.80 X 10^ 


6.76 X 10^ 




I'raA 


cm^ s 


1.74 X 10^ 


1.43 X 10^ 




V 


2 —1 

cm s 


19.3 


6.90 X 10^ 




K 


2 — 1 

cm g 


0.382 


0.728 




\f\ 


cms"^ 


3.91 X 10-5 


8.50 X 10-5 




n 


rads"^ 


10-3 


lo-'^ 





Table 2. Magnetic Flux Ring Parameters 



Our Model (A/x = 5 x IQ-^) BWNC Model (A/x = 0) 



Parameter 


Units 


Bf = 


7.3 G 


Bj- 


= 73G 


RGB-1 


RGB-2 


ao 


cm 


1.3 X 


10^ 


9.7 X 


10*^ 


6.5 X 10*^ 


1.5 X 10^ 


ao 


Hp 


1.5 X 


10-2 


1.1 X 


10-2 


5.1 X 10-3 


1.2 X 10-2 


tb 


yr 


2.9 X 


102 


1.6 X 


102 


0.19 


1.1 






7.3 




13 




9.8 X 10^ 


1.8 X 10^ 


(-B^)o 


kG 


254 




254 




380 


48 


/3o 




1.8 X 


10^ 


1.8 X 


10^ 


3.3 X 10*^ 


2.0 X 10^ 


(■B<^)bce 


kG 


3.5 




3.5 




3.5 


0.44 


N 




1.3 X 


102 


1.3 X 


102 


1.0 


1.0 


N = N/tb 


yr-i 


0.45 




0.83 




5.3 


0.91 



